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Cosmological observations in a modified theory of gravity (MOG) 
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Our modified gravity theory (MOG) is a gravitational theory without exotic dark matter, 
based on an action principle. MOG has been used successfully to model astrophysical 
phenomena such as galaxy rotation curves, galaxy cluster masses, and lensing. MOG may 
also be able to account for cosmological observations. We assume that the MOG point 
source solution can be used to describe extended distributions of matter via an appro- 
priately modified Poisson equation. We use this result to model perturbation growth in 
MOG and find that it agrees well with the observed matter power spectrum at present. 
As the resolution of the power spectrum improves with increasing survey size, however, 
significant differences emerge between the predictions of MOG and the standard ACDM 
model, as in the absence of exotic dark matter, oscillations of the power spectrum in 
MOG are not suppressed. We can also use MOG to model the acoustic power spectrum 
of the cosmic microwave background. A suitably adapted semi-analytical model offers 
a first indication that MOG may pass this test, and correctly model the peak of the 
acoustic spectrum. 

Keywords: Cosmology:theory — Large-scale structure of Universe — Gravitation. 



1. Introduction 

The preferred model of cosmology today, the ACDM model, provides an excellent 
fit to cosmological observations, but at a substantial cost: according to this model, 
about 95% of the universe is either invisible or undetectable, or possibly boM. This 
fact provides a strong incentive to seek alternative explanations that can account 
for cosmological observations without resorting to dark matter or Einstein's cosmo- 
logical constant. 

For gravitational theories designed to challenge the ACDM model, the bar is set 
increasingly higher by recent discoveries. Not only do such theories have to explain 
successfully the velocity dispersions, rotational curves, and gravitational lensing of 
galaxies and galaxy clusters, the theories must also be in accord with cosmolog- 
ical observations, notably the acoustic power spectrum of the cosmic microwave 
background (CMB), the matter power spectrum of galaxies, and the recent obser- 
vation of the luminosity-distance relationship of high-z supernovae, which is seen 
as evidence for "dark energy" . 

Modified Gravity (MOG 2 ) has been used suc cessfully to account for galaxy 
cluster masse the rotation curves of galaxieJ^^l, velocity dispersions of satellite 
galaxies^, andglobular clusters^. It was also used to offer an explanation for the 
Bullet Cluster without resorting to nonbaryonic dark matter. 

MOG may also be able to meet the challenge posed by cosmological observations. 
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We investigate two sets of observations in particular: the matter power spectrum 
that describes the spatial distribution of galaxies in the Universe, and the acoustic 
spectrum of the cosmic microwave background (CMB) radiation. 

In the next section, we review the key features of MOG. This is followed by 
sections presenting detailed calculations for the galaxy power spectrum and the 
acoustic power spectrum of the CMB. A concluding section summarizes our results 
and maps out future steps. 



2. Modified Gravity Theory 

Modified Gravity (MOG) is a fully relativistic theory of gravitation that is derived 

191 

from a relativistic action principle involving scalar, tensor, and vector fields. MOG 
has evolved as a result of investigations of Nonsymmetric Gravity Theory (NGT^J), 
and most recently, it has taken the form of Scalar- Tensor- Vector Gravity (STVG=|). 
In the weak field approximation, STVG, NGT, and Metric-Skew- Tensor Gravity 
(MSTGpl) produce similar results. 



2.1. Scalar-Tensor-Vector Gravity 



Our modified gravity theory is based on postulating the existence of a massive vector 



field, 



The choice of a massive vector field is motivated by our desire to introduce 



a repulsive modification of the law of gravitation at short range. The vector field is 
coupled universally to matter. The theory, therefore, has three constants: in addition 
to the gravitational constant G, we must also consider the coupling constant u that 
determines the coupling strength between the field and matter, and a further 
constant fi that arises as a result of considering a vector field of non-zero mass, and 
controls the coupling range. The theory promotes G, /z, and lo to scalar fields, hence 
they are allowed to run, resulting in the following actiorJ^JJ 



S — Sg + Si/, + Ss + S 



m • 
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(3) 

(4) 



Here, 

Vju 



Sm is the "matter" action, while = d^4> y — d v <\>^ and ((f)), Vg(G), 



and Vjj,(fi) denote the self- interaction potentials associated with the vector 
field and the three scalar fields. The symbol V M is used to denote covariant differen- 
tiation with respect to the metric g^ v , while the symbols R, A, and g represent the 
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Ricci-scalar, the cosmological constant, and the determinant of the metric tensor, 
respectively. We define the Ricci tensor as 

R^v = d a v °„ — d v Y " Q + r™„r aj a — r^rf „. (5) 

Our units are such that the speed of light, c = 1; we use the metric signature 
(+,-,-,-)• 

The apparent "wrong" sign of the V^wVw term in the Lagrangian is of po- 
tential concern; however, we found that in all the solutions (including numerical 
solutions) considered to date, u remains constant. Keeping w as a dynamical scalar 
field (with the "wrong" sign in the Lagrangian) allowed us to develop a parameter- 
free solutiorlHl, but we anticipate that the ui field may disappear from the theory 
as it is being further developed. 

A direct numerical solution of the theory's field equations in the spatially homo- 
geneous, isotropic case (FLRW cosmology) yields an expanding universe. Choosing 
a constant Vq as one of the initial parameters of the solution, the age of the universe 
can be adjusted to fit observation. As an alternative, we also considered changing 
the overall sign of the kinetic terms in Ss', this solution, which violates several en- 
ergy conditions, but keeps the energy density p positive, is a "bouncing" cosmology 
(indeed, a classical bouncing cosmology requires that some or all of the energy con- 
ditions be violated) . In this cosmology, the age of the universe since the bounce and 
the density of the universe at the time of the bounce can be tuned by choosing an 
appropriate constant Vq. Either way, a solution in which the age of the universe 
is in agreement with observation can be obtained. These solutions are a subject of 
further study, which will be reported elsewhere. 



2.2. Point particles in a spherically symmetric field 

For a point particle moving in the spherically symmetric field of a gravitating source, 
a particularly simple solution for the acceleration is obtained^. 

r = -^T- [l + a-a(l + He-" r ] , (6) 

where M is the source mass, while a determines the strength of the "fifth force" 
interaction, and fi controls its range. In prior work, a and /i were considered free 
parameters that were fitted to data. Our recent wor allows us to determine a 
and /i as functions of the source mass M: 

M /G "-lV (7) 



(VM + £) 2 \G n 
and 



This solution can be seen to satisfy the field equations in the spherically symmet- 
ric case either numerically or by deriving an approximate solution analytically^! 
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The numerical values for D and E are determined by matching the result against 
galaxy rotation curves^} 

D ~ 6250 M^kpc -1 , (9) 

E ~ 25000 A4 /2 - (10) 

The value of Goo ~ 20Gn is set to ensure that at the horizon distance, the effective 
strength of gravity is about 6 times Gjv, eliminating the need for cold dark matter 
in cosmological calculations, as described in the previous section. 



2.3. The MOG Poisson Equation 

The acceleration law ([6]) is associated with the potential, 

a 



$ = _G^M 



1 - 



l + a 



= $ N + $ Y , (11) 



where 



$at = (12) 

r 

is the Newtonian gravitational potential with Goo = (1 + a)Gjy as the gravitational 
constant, and 



$y = -—G^M (13) 

l + a r 

is the Yukawa-potential. These potentials are associated with the corresponding 
Poisson and inhomogeneous Helmholtz equations, which are given by Ref. |5J 

V 2 $jv(r) = 4 7 rG 00/9 (r), (14) 

(V 2 - M 2 )<Mr) - -4^-^G ooP (r). (15) 
l + a 

Full solutions to these potentials are given by 

$jv(r) = -Goo / d 3 r, (16) 
J |r-r| 

*r W - ^Goo / d^r. (17) 

l + a J \r — t\ 

These solutions can be verified against Eqs. (p~2]) and ([T3l) by applying the delta 
function point source density p(r) = M<5 3 (r). 

Strictly speaking, (|17[) is a valid solution only when a is approximately constant. 
For inhomogeneous matter distributions, a is expected to vary as a function of 
matter density, and as such, this naive application of the spherically symmetric, 
static vacuum solution to model extended distributions of matter breaks down. 
However, for small perturbations of a homogeneous background, we expect (|17j) 
to remain valid; this expectation can be verified once generalized (approximate) 
solutions of the theory in the presence of matter become available. 
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Combining Eq. (JTTJ) with Eqs. flU} and (JI5| yields 
V 2 $ = 47rGjv/3(r) + p. 2 ^r{r) 

= 47rGAr/9(r) + a^ 2 G N / : — 




(18) 



containing, in addition to the usual Newtonian term, a nonlocal source term on the 
right-hand side. 

3. MOG and the matter power spectrum 

The distribution of mass in the universe is not uniform. Due to gravitational self- 
attraction, matter tends to "clump" into ever denser concentrations, leaving large 
voids in between. In the early universe, this process is counteracted by pressure. 
The process is further complicated by the fact that in the early universe, the energy 
density of radiation was comparable to that of matter. 

3.1. Density fluctuations in Newtonian gravity 

To first order, this process can be investigated using perturbation theory. Taking 
an arbitrary initial distribution, one can proceed to introduce small perturbations 
in the density, velocity, and acceleration fields. These lead to a second-order dif- 
ferential equation for the density perturbation that can be solved analytically or 
numerically. This yields the transfer function, which determines how an initial den- 
sity distribution evolves as a function of time in the presence of small perturbations. 

3.1.1. Newtonian theory of small fluctuations 

In order to see how this theory can be developed for MOG, we must first review how 
the density perturbation equation is derived in the Newtonian case. Our treatment 
follows closely the approach presented by Ref. [13] We begin with three equations: 
the continuity equation, the Euler equation, and the Poisson equation. 




0. 



4nGp. 



-Vp + g, 



P 



(19a) 



(196) 



(19c) 



First, we perturb p, p, v and g. Spelled out in full, we get: 




+ [(v + 5v) • V](v + 5v) = ^V(p + 5p) + g + 6g, 

p + 5p 

V ■(s + 5g) = -ATrG{p + 5p). 




(20a) 



(206) 



(20c) 
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Subtracting the original set of equations from the new set, using l/(p + Sp) = 
(p — 8p)/[p 2 — (Sp) 2 ] = 1/p— Sp/ p 2 , and eliminating second-order terms, we obtain 

^+V-(<5pv + p<5v) = 0, (2l a ) 

— ^ + (v V)ov + (5v • V)v = -?Vp - -VSp + Sg, (216) 
ot p z p 

V • 5g = -AirGSp. (21c) 

A further substitution can be made by observing that Sp = (Sp/Sp)Sp = c 2 s 8p 
where c 2 — (dp /dp) adiabatic is the speed of sound. We can also eliminate terms 
by observing that the original (unperturbed) state is spatially homogeneous, hence 
Vp = Vp = 0: 

^+v- VSp + 5pV-v + pV-Sv = 0, (22a) 
BSw c 2 

— + (v V)<fv + (5v • V)v = --^VSp + 5g, (22b) 
V • Sg = -4wGSp. (22c) 

Now we note that v = Hx, hence 

V • v = HV • x = 3H, 
(Sv ■ V)v = (<5v • V)(iJx) = H(Sw ■ V)x = HSw. 

Therefore, 

^ + v VSp + SHSp + pV ■ 6v = 0, (23a) 
BSv c 2 

— + (v-V)Sv + HSv=-^VSp + Sg, (23b) 
Ot p 

V • 5g = -4nG5p. (23c) 

The next step is a change of spatial coordinates to coordinates comoving with the 
Hubble flow: 

x = a(t)q. 

This means 

d\ (d_ 
Ml „ ~ \dt 

and 



+ vV x , 



V q = aV x . 
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After this change of coordinates, our system of equations becomes 

^ + 3H6p+-pV -5v = 0, (27a) 
at a 

BS~v c 2 

— +H5v = — s -VSp + Sg, (276) 
at ap 

V • 5g = -AnaGSp. (27c) 

Now is the time to introduce the fractional amplitude S = Sp/p. Dividing (I27a|) 
with p, we get 

5+^S + 3HS+-V-Sv = Q. (28) 
p a 

However, since p = pod^/a 3 , and hence p/p = — 3d/a, the second and third terms 
cancel, to give 

-aS = VSv. (29) 
Taking the gradient of (|276[) and using (|27c| to express V • <5g, we get 



8 r 2 

— (-aS) +H(-aS) = --^V 2 5-^Ga P 5. (30) 

Spelling out the derivatives, and dividing both sides with a, we obtain 

r 2 

S + 2HS - -|V 2 £ - AirGpS = 0. (31) 
a z 

For every Fourier mode 6 = S^(t)e lk ci (such that V 2 <5 = ~k 2 S), this gives 

J k + 2i?(5 k + f - AnGp 4 = 0. (32) 

The quantity fc/a is called the co- moving wave number. 

If fc is large, solutions to ([3"2l are dominated by an oscillatory term; for small k, 
a growth term predominates. 

A solution to (|3"2"]) tells us how a power spectrum evolves over time, as a function 
of the wave number; it docs not specify the initial power spectrum. For this reason, 
solutions to (I32[) are typically written in the form of a transfer function 

T(k) = ^ = 0)^-oo) 
6 k {z = oo)d {z = 0) 

If the initial power spectrum and the transfer function are known, the power 
spectrum at a later time can be calculated (without accounting for small effects) as 

P(k)=T 2 (k)P (k). (34) 

P(k) is a dimensioned quantity. It is possible to form the dimensionless power 
spectrum 

A 2 (k) = Ak 3 T 2 (k)P (k), (35) 
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where A is a normalization constant determined by observation. This form often 
appears in the literature. In the present work, however, we are using P(k) instead 
of A(Jfe). 

The initial power spectrum is believed to be a scale invariant power spectrum: 

P (k) oc k n , (36) 
where n ~ 1. A recent estimation n is n = 0.963lg'Q}g. 



3.1.2. Analytical approximation 



Eq. (|32j) is not difficult to solve in principle. The solution can be written as the 
sum of oscillatory and growing terms. The usual physical interpretation is that 
when pressure is sufficient to counteract gravitational attraction, this mechanism 
prevents the growth of density fluctuations, and their energy is dissipated instead 
in the form of sound waves. When the pressure is low, however, the growth term 
dominates and fluctuations grow. Put into the context of an expanding universe, 
one can conclude that in the early stages, when the universe was hot and dense, the 
oscillatory term had to dominate. Later, the growth term took over, the perturbation 
spectrum "froze" , affected only by uniform growth afterwards. 

In practice, several issues complicate the problem. First, the early universe can- 
not be modeled by matter alone; it contained a mix of matter and radiation (and, 
possibly, neutrinos and cold dark matter.) To correctly describe this case even using 
the linear perturbation theory outlined in the previous sections, one needs to re- 
sort to a system of coupled differential equations describing the different mediums. 
Second, if the perturbations are sufficiently strong, linear t heo ry may no longer be 
valid. Third, other nonlinear effects, including Silk-damping^^, cannot be excluded 
as their contribution is significant (indeed, Silk damping at higher wave numbers is 
one of the reasons why a baryon-only cosmological model based on Einstein's theory 
of gravity without dark matter fails to account for the matter power spectrum.) 

The authors of Ref. [15] addressed all these issues when they developed a semi- 
analytical solution to the baryon transfer function. This solution reportedly yields 
good results in the full range of < fib < 1. Furthermore, unlike other approx- 
imations and numerical software codes, this approach keeps the essential physics 
transparent, allowing us to adapt the formulation to the MOG case. 

In Ref. [15] the transfer function is written as the sum of a baryonic term Tb and 
a cold dark matter term T c : 

T(k) = ^T b (k) + ^T c (k), (37) 

where fi c represents the cold dark matter content of the universe relative to the crit- 
ical density. As we are investigating a cosmology with no cold dark matter, we ignore 
T c . The baryonic part of the transfer function departs from the cold dark matter 
case on scales comparable to, or smaller than, the sound horizon. Consequently, the 
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baryonic transfer function is written as 

fo(M.l) , "be- 



T b (k) 



with 



where 



and 



1 + (ks/5.2) 2 1 + {pb/ksf 
In (e + 1.8/3 c ?) 



sin ks 
ks 



T (fc,a c ,/3 C ) 



C 



In (e + 1.8/3 c g) + Cq 2 ' 

14.2 386 
~ + l + 69.9g 108 ' 



q = k@l 7 (n m h 2 ) \ 

The sound horizon is calculated as 



_ 2 / G ^ V 1 + -Rd + \jRd + Rcq 
3fceq y Rcq 



1 + \/Req 

The scale at the equalization epoch is calculated as 

fc cq = 7.46 x 10" 2 O m /i 2 97T. 



'm' 1 ^2.7- 



(38) 

(39) 
(40) 
(41) 

(42) 
(43) 



The transition from a radiation-dominated to a matter-dominated era happens at 
the redshift 



z cq = 25000il m h 2 Q^, 
while the drag era is defined as 

(n ro /^ - 251 



where 



and 



^ = 1291 l + 0.65V m ^ 8 [1 + 6l(nm;t2)ta 



h = 0.313(Q m /i 2 )- U419 [l + 0.607(Q m /i 2 ) Ub74 ], 



b 2 = 0.238(fi m /i 2 ) - 223 . 
The baryon-to-photon density ratio at a given redshift is calculated as 



R = 31.5fi m /i 2 9^7 



1000 



(44) 

(45) 

(46) 
(47) 

(48) 



The Silk damping scale is obtained using 

fcsiik = l.6(nbh 2 )°- 52 (n m h 2 ) a - 73 [l + (10.4Q™/* 2 )" - 95 ]. (49) 
The coefficients in the second term of the baryonic transfer function are written as 

2.07fc cqS (l + i? d )- 3 / 4 F(i±^V (50) 

V 1 + z d J 
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Fig. 1. The matter power spectrum. Three models are compared against five data sets (see text): 
ACDM (dashed blue line, Sl b = 0.035, Q c = 0.245, Q A = 0.72, H = 71 km/s/Mpc), a baryon-only 
model (dotted green line, Q b = 0.035, H = 71 km/s/Mpc), and MOG (solid red line, a = 19, 
H = 5h Mpc" 1 , U b = 0.035, H = 71 km/s/Mpc.) Data points are colored light blue (SDSS 2006), 
gold (SDSS 2004), pink (2dF), light green (UKST), and dark blue (CfA). 



fa = 0.5 + -± + 3 - 2-^ ) y/(17.2Sl m h*)2 + 1, 



where we used the function 



F(y) = y 



-6 v /TT^+(2 + 3 2 /)ln 



1 



(51) 



(52) 



A shifting of nodes in the baryonic transfer function is accounted for by the quantity 



m = 



1/3' 



(53) 



1 + (Pnodc/ks)" 

where 

(3 nodc = 8Al(n m h 2 ) 0A35 . (54) 

The symbol Q 2 .7 = T/2.7 is the temperature of the CMB relative to 2.7 K, while 
h = iJ/(100 km/s/Mpc). The wave number k is in units of Mpc -1 . 



3.2. Density fluctuations in Modified Gravity 

We consider the MOG Poisson equation (jT5|) . established in section 12.31 As the 
initial unperturbed distribution is assumed to be homogeneous, p is not a function 
of r and can be taken outside the integral sign: 



(55) 



19, 2013 20:51 WSPC/INSTRUCTION FILE ms 



11 



Varying p, we get 



V • 5g(r) = -AirG N 5p{v) 



-p Gn&5p 



r - r' 



- M |r-r'| d 3 r /_ 



Accordingly, (T5TI) now reads 



S + 2H5 - -|V 2 <5 - AttG nP S - p 2 G N apS 



,-/i|r-r 



-rfV = 0. 



The integral can be readily calculated. Assuming that |r — r'| runs from 
comoving wavelength a/k, we get 



(56) 

(57) 
to the 



,-/i|r-r 



tt/2 2ir a/k 



3„' 



-cfr 



r sin 6* dr dd> d6 



ooo 
4tt [l - (1 + pa/k)e' tia,k } 
7 



Substituting into (1571) . we get 



S + 2H5 - -^V 2 5 - 4nG N pS 



AttGnch 



1 - ! 1 



P s = o, 



or 



5 + 2H6--±V 2 6 



-4nG N {l + a [l - (l + ^) e -Wfc] J pg = o. 

This demonstrates how the effective gravitational constant 
G ea = G iV {l + a[l-(l + ^)e-^]} 

depends on the wave number. 

Using G e ff, we can express the perturbation equation as 



iirGcsp <5 k = 0. 



<5 k + 2£T<5 k - 

As the wave number k appears only in the source term 

'c 2 k 2 



(58) 



(59) 
(60) 



(61) 



(62) 



47rG c ffP , 



it is easy to see that any solution of (|3"2"]) is also a solution of (|6"2"j) . provided that k 
is replaced by k! in accordance with the following prescription: 



, , 2 //2 , ( /(? off- g AT > A _ 2 

Gat 



(64) 
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where Aj = y 'c^/GnP is the Jeans wavelength. 

This shifting of the wave number applies to the growth term of the baryonic 
transfer function (|38|) . However, as the sound horizon scale is not affected by changes 
in the effective gravitational constant, terms containing the product ks must remain 
unchanged. Furthermore, the Silk damping scale must also change as a result of 
changing gravity; this change is proportional to the 3/4th power of G, as demon- 
strated bjfl Ref. HU thus 

/Q cff \ 3/4 

fc snk = fc siik f -^-J , (65) 

(note also Eq. H9"|) . Using these considerations, we obtain the modified baryonic 
transfer function 

f (fc',l,l) a b cxp(-[fc/^ ilk ] 1 - 4 ) | 
1 + (Ls/5.2) 2 + 1 + {pb/ksf J ' K ' 

The effects of these changes can be summed up as follows. At low values of k, 
the transfer function is suppressed. At high values of k, where the transfer function 
is usually suppressed by Silk damping, the effect of this suppression is reduced. The 
combined result is that the tilt of the transfer function changes, such that its peaks 
are now approximately in agreement with data points, as seen in Figure [TJ 

Data points shown in this figure come from several source s. Firs t and foremost, 
the two data releases of the Sloan Digital Sky Survey 

(SDSS^im 

are presented. 

Additionally, data from the 2dF Galaxy Redshift SurvejffS, UKSlP \ and CfA13(Pl 
surveys are shown. Apart from normalization issues, the data from these surveys are 
consistent in the range of 0.01 h Mpc -1 < k < 0.5 h Mpc -1 . Some surveys provide 
data points outside this range, but they are not in agreement with each other. 

3.3. Discussion 

As a result of the combined effects of dampened structure growth at low values of 
k and reduced Silk damping at high values of k, the slope of the MOG transfer 
function differs significantly from the slope of the baryonic transfer function, and 
matches closely the observed values of the matter power spectrum. On the other 
hand, the predictions of MOG and ACDM cosmology differ in fundamental ways. 

First, MOG predicts oscillations in the power spectrum, which are not smoothed 
out by dark matter. These oscillations may be detectable in future galaxy surveys 
that utilize a large enough number of galaxies, and sufficiently narrow window 
functions in order to be sensitive to such fluctuations. However, the finite size of 
samples and the associated window functions used to produce presently available 
power spectra mask any such oscillations. To illustrate this, we applied the same 
window function to the MOG prediction, which resulted in a smoothed curve, seen in 



a cf. Eq. (4.210) in Rel.E] note that fi/i 2 <x G. 
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P(k) 




5 X 10 3 



Fig. 2. The effect of window functions on the power spectrum is demonstrated by applying the 
SDSS luminous red galaxy survey window functions to the MOG prediction. Baryonic oscillations 
are greatly dampened in the resulting curve (solid red line). A normalized linear ACDM estimate 
is also shown (thin blue line) for comparison. 

Figure [2] A \ 2 comparison actually suggests that MOG offers a better fit (xmog = 
0.03, 

Xacdm = 0-09 P er degrees of freedom), although we must be cautious: the 
ACDM approximation we used is not necessarily the best approximation available, 
and the MOG result is dependent on the validity of the analysis presented in this 
section, which was developed without the benefit of an interior solution. 

Second, MOG predicts a dampened power spectrum at both high and low values 
of k relative to ACDM. Observations at sufficiently high values of k may not be 
practical, as we are entering sub-galactic length scales. Low values of k are a different 
matter: as accurate three-dimensional information becomes available on ever more 
distant galaxies, power spectrum observations are likely to be extended in this 
direction. 

In the present work, we made no attempt to account for the possibility of a non- 
zero neutrino mass, and its effects on the power spectrum. Given the uncertainties 
in the semi-analytical approximations that we utilized, such an attempt would not 
have been very fruitful. Future numerical work, however, must take into account 
the possibility of a non-negligible contribution of neutrinos to the matter density. 



4. MOG and the CMB 

The cosmic microwave background (CMB) is highly isotropic, showing only small 
temperature fluctuations as a function of sky direction. These fluctuations are not 
uniformly random; they show a distinct dependence on angula r si ze, as has been 
demonstrated by the measurements of the Boomerang experiment!^ and the Wilkin- 
son Microwave Anisotropy Probe 

The angular power spectrum of the CMB can be calculated in a variety of ways. 
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The preferred method is to use numerical software, such as CMBFASlP^l. Unfortu- 
nately, such software packages cannot easily be adapted for use with MOG. Instead, 
at the present time we opt to use the excellent semi-analytical approximation de- 
veloped by Ref. 23, While not as accurate as numerical software, it lends itself more 
easily to nontrivial modifications, as the physics remains evident in the equations. 

What justifies the use of this semi-analytical approach is the fact that the phe- 
nomenology of MOG vs. dark matter can be understood easily. Collisionless cold 
dark matter interacts with normal matter only through gravity. In the late universe, 
the ratio of cold dark matter vs. baryonic matter varies significantly from region to 
region; this is why the results of the previous section are nontrivial and significant. 
However, in the early universe (recombination era), the universe was still largely 
homogeneous, and cold dark matter effectively acted as a "gravity enhancer": its 
effects can be mimicked by simply increasing the effective gravitational constant. 

This may seem surprising in view of studies that have placed stringent con- 
straints on the variability of G. For example, after substituting G — > X 2 G in the 
Friedmann equation, the authors of Ref. I2"4"lhave shown that A is constrained to 10% 
or better by WMAP data. At first sight, this seems inconsistent with our assertion 
that MOG, with G c s > Gn, can successfully mimic the effects of dark matter on 
the CMB acoustic spectrum. Yet that this is the case can be seen if one writes down 
the Friedmann equation after incorporating A: 



The full form of the substitution rule, therefore, is Gp — > X 2 Gp. In MOG, we sub- 
stitute G — > G c ff and p — > pi, (no CDM component), but (Gp) acdm = (G c gpb)MOGj 
hence A = 1. 

This discussion leads to a simple substitution rule that is applicable when the 
universe is approximately homogeneous. When a quantity containing G appears 
in an equation describing a gravitational interaction, G c ff must be used. However, 
when a quantity like fij, is used to describe a nongravitational effect, the Newtonian 
value of Gn must be retained. 

Our choice to use Mukhanov's semianalytical approximation is motivated by the 
fact that these substitutions can be made in the formulae in a straightforward and 
unambiguous manner. 

4.1. Semi- analytical estimation of CMB anisotropics 

In Ref. [23] we find a calculation of the correlation function C(l), where I is the 
multipole number, of the acoustic power spectrum of the CMB using the solution 



(67) 
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where I ^ 1,0 denotes the oscillating part of the spectrum, while the non-oscillating 
part is written as the sum of three parts: 



N = N x + N 2 + N 3 . 
These, in turn, are expressed as 



JVi = 0.063e 
N 2 



2 [P - 0.22(1/1 f ) - 6 -2.6] 
1 + 0.65(1/1 f ) 1A 

0.037 [P-0.22(V^) a3 + 1.7] 2 



(1 + £)V2 l + 0.65(l/l s ) 1A 
0.033 [P-0.5(V^)°' 55 + 2.2] 2 



The oscillating part of the spectrum is written as 
O = 



-d/iff 



(68) 

(69) 
(70) 
(71) 



A\ cos (pi + ^ + A 2 cos (2 pi 



where 



and 



A 1 = 0.1£ 



(P-0.78) 2 -4.3 1 v 



e 2 



A 2 = 0.14 



(1 + 1/4 

(0.5 + 0.36P) 2 



(72) 
(73) 
(74) 



(1+0 1/2 ' 

The parameters that occur in these expressions are as follows. First, the baryon 
density parameter: 

£ = 17 (n b h 2 75 ) , (75) 

where ~ 0.035 is the baryon content of the universe at present relative to the 
critical density, and /175 = H/(75 km/s/Mpc). The growth term of the transfer 
function is represented by 

C7 0.09/ 

P = In m -J =, ( 76 ) 

200^/0^4 

where Sl m ~ 0.3 is the total matter content (baryonic matter, neutrinos, and cold 
dark matter). The free-streaming and Silk damping scales are determined, respec- 
tively, by 



If = 1300 [l + 7.8 x 10~ 2 (n m h 2 5 ) 1 



1/2 



0.09 



0.7/, 



l+0.56g 



0.8 



("->js)' 



(77) 
(78) 



1+ « [i+(i+^n m ^ 5 )- 1 ^]- 
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Fig. 3. MOG and the acoustic power spectrum. Calculated using !!m = 0.3, Qb = 0.035, Ho = 
71 km/s/Mpc. Also shown are the raw WMAP 3-year data set (light blue), binned averages 
with horizontal and vertical error bars provided by the WMAP project (red), and data from the 
Boomerang experiment (green). 



Lastly, the location of the acoustic peaks is determined by the parameteiQ 

p = 0.015(1 + O.UO^ittmhfs 1 ) 016 . (79) 



4.2. The MOG CMB spectrum 

The semi-analytical approximation presented in the previous section can be adapted 
to the MOG case by making two important observations. 

First, in all expressions involving the value of Mukhanov's f2 m (which includes 
contributions from baryonic matter and cold dark matter using Newton's gravita- 
tional constant), we need to use — 0.3 (which includes baryonic matter only, 
using the running value of the gravitational constant, G e ff — GGn)- Second, we 
notice that the value of Of, in (|75p does not depend on the effective value of the 
gravitational constant, as this value is a function of the speed of sound, which de- 
pends on the (baryonic) matter density, regardless of gravitation. In other words, 
fl b ~ 0.035 is calculated using Newton's gravitational constant. 

After we modify Mukhanov's semi- analytical formulation by taking these con- 
siderations into account, we obtain the fit to the acoustic power spectrum shown in 
Figure [3] 



b Note that we slightly adjusted the coefficients of 177> and J79j, which improved the fit noticeably, 
while remaining fully consistent with Mukhanov's derivation. 
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4.3. Discussion 

As Figure [3] demonstrates, to the extent that Mukhanov's formulation is applicable 
to MOG, the theory achieves agreement with the observed acoustic power spectrum. 
This result was obtained without fine-tuning or parameter fitting. The MOG con- 
stant fj, was assumed to be equal to the inverse of the radius of the visible universe. 
Thereafter, the value of a is fixed if we wish to ensure — 0.3. This was sufficient 
to achieve consistency with the data. 



5. Conclusions 

In this paper, we demonstrated how MOG can account for key cosmological obser- 
vations using a minimum number of free parameters. We applied the MOG point 
source solution in a suitably modified form of the Poisson equation and re-derived 
the equations of structure growth. We found that the result is in agreement with 
presently available observational data. 

Notably, we also found that as the available data sets grow in size, a significant, 
and likely irreconcilable disagreement emerges between the predictions of MOG 
and those of the ACDM concordance model. In ACDM, the presence of collisionless 
exotic dark matter leads to a significant dampening of the baryonic oscillations in 
the matter power spectrum: unit oscillations are suppressed, and appear only as a 
slight modulation of the power spectrum at shorter wavelengths. In contrast, unit 
oscillations are not suppressed in MOG. Presently, these oscillations are not seen 
only because the resolution of the data is not high enough: when we apply the 
appropriate bin sizes and window functions to a simulated data set, the resulting 
curve is nearly smooth. As galaxy surveys grow in size, however, bin sizes will get 
smaller and, if MOG is correct, the unit oscillations will emerge in the data. 

We also investigated the acoustic power spectrum of the cosmic microwave back- 

IOOI 

ground using MOG. Existing software codes, notably the program CMBFAST and 
its derivatives, are ill suited for this investigation as it is difficult to disentangle 
the use of quantities proportional to Gp in gravitational vs. nongravitational con- 
texts. Before embarking on what seems to be a formidable task, we turned to a 
semi- analytical approximation^! While many of the approximations employed by 
Ref. [23] are not physically motivated, but numerical fitting formulae, nonetheless 
the role played by quantities proportional to Gp can be clearly discerned, and the 
formulae can be suitably adapted. While we recognize that this is not a conclusive 
result, we find it nonetheless encouraging that the CMB acoustic power spectrum 
was faithfully reproduced. 

In conclusion, we have demonstrated that cosmological observations of the mat- 
ter power spectrum and the CMB acoustic spectrum do not trivially rule out MOG 
as a possible alternative to the standard ACDM model of cosmology. 
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